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Abstract 

Phase reddening is an effect that produces an increase of the spectrai slope and variations 
in the strength of the absorption bands as the phase angle increases. In order to understand 
its effect on spectroscopic observations of asteroids, we have analyzed the visible and near- 
infrared spectra (0.45-2.5 /jm) of 12 near-Earth asteroids observed at different phase angles. 
All these asteroids are classified as either S-complex or Q-type asteroids. In addition, we 
have acquired laboratory spectra of three different types of ordinary chondrites at phase 
angles ranging from 13° to 120°. We have found that both, asteroid and meteorite spectra 
show an increase in band depths with increasing phase angle. In the case of the asteroids 
the Band I depth increases in the range of ~ 2° < g < 70° and the Band II depth increases 
in the range of ~ 2° < g < 55°. Using this information we have derived equations that can 
be used to correct the effect of phase reddening in the band depths. Of the three meteorite 
samples, the (olivine-rich) LL6 ordinary chondrite is the most affected by phase reddening. 
The studied ordinary chondrites have their maximum spectral contrast of Band I depths at a 
phase angle of ~ 60°, followed by a decrease between 60° and 120° phase angle. The Band II 
depths of these samples have their maximum spectral contrast at phase angles of 30° — 60° 
which then gradually decreases to 120° phase angle. The spectral slope of the ordinary 
chondrites spectra shows a significant increase with increasing phase angle for g > 30°. 
Variations in band centers and band area ratio (BAR) values were also found, however they 
seems to have no significant impact on the mineralogical analysis. Our study showed that 
the increase in spectral slope caused by phase reddening is comparable to certain degree of 
space weathering. In particular, an increase in phase angle in the range of 30° to 120° will 
produce a reddening of the reflectance spectra equivalent to exposure times of ~ 0.1 x 10 6 
to 1.3 x 10 6 years at about 1 AU from the Sun. This increase in spectral slope due to phase 
reddening is also comparable to the effects caused by the addition of different fractions of 
SMFe. Furthermore, we found that under some circumstances phase reddening could lead 
to an ambiguous taxonomic classification of asteroids. 



*E-mail address: sanchez@mps.mpg.de 
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1 Introduction 



Near-infrared (NIR) spectroscopy is a powerful technique to derive information about the surface 
mineralogy of asteroids. This mineralogical characterization relies primarily on the identifica- 
tion and analysis of diagnostic features present in the spectra of some asteroids. Olivine and 
pyroxene are two common mafic minerals found on asteroids ( Gaffev et al. . 20021 ) . and their 



spectral properties dominate the reflect ance spectra of som e particular classes. Asteroids be- 
longing to the S-complex and Q- types ( DeMeo et al. . 20091 . Bus-DeMeo taxonomy hereafter) 
are examples of these objects. In the case of olivine, the absorption feature is composed of three 
overlapping bands and is centered near 1.04-1.1 //m, while pyroxenes show to major absorp- 
tion bands, c entered near 0.9-1 um and 1.9-2 /im, all of them caused by the presence of Fe 2+ 
cations (e.g., Adams . 19741 . 19751 : Burns . 1993 ). In the spectra of olivine- pyroxene mixtures the 
wavelength position (band centers) of the combined absorption features near 1 /xm (Band I) is 
a function of relative abundance and composition of olivine and pyroxene, while the posi t ion of 



the feature near 2 / /m (Band II) is a function of the pyroxene composition ( Cloutis et al. . 19861 ; 



Gaffev et allll993h . In addition to the absorption band centers, the ratio of the area of Band 



II to that of Band I, known as band area ratio (BAR), is used to estimate o l ivine and pyroxene 



abun dances in asteroids and meteorites (jCloutis et all Il986l ; iBurbine et all 120031 : iDunn et al 



2010h . Together, these spectral band parameters constitute a useful tool to obtain information 
about the surface mineralogy and composition of asteroids. However, band parameters can 
also be affected by other factors not related to compositional variations. One of these factors 
is the phase angle at the moment of the observation. The phase angle (g), is defined as the 
angular separation between the Sun and the observer as seen from the target. Phase angle in- 
duced effects can manifest themselves as phase reddening, which produces an artificial increase 
(reddening) of the spectral slope and variations in the strength of the absorption bands with 
increasing phase angle. Thi s effect is explained as the result of the wavelength dependence o f 
the single-scattering albedo ( Gradie et al. . 1980l ; Gradie and Veverka . 1986 : Clark et al. . 2002bl ). 
Traditionally, the term phase reddening was used to describe only the increase of the spectral 
slope or continuum, however as will be shown later, as the phase angle increases variations in 
the absorption bands are also seen. Therefore, in the present work the term phase reddening is 
extended to refer not only to the increase of the spectral slope, but also to the variations in the 
strength of the absorption bands. 

Ph ase reddening e ffect among Solar System objects was first noticed in asteroid broad band 
colors. Gehrels ( 1970l ) reported a phase reddening in the B-V and U-B colors of 4 Vesta to be 
0.0018 and 0.0027 m ag/degree, re s pectiv ely. Additional studies us i ng ph otometric observations 
include the works of Millis et al. ( 19761 ) and Lumme and Bowelll dl98lh . The effects of phase 
reddening have been also observed among NEAs. iLuu and Jewittl (|l990l ) found that the spectral 
slopes of NEAs in general are higher than those measured for 3:1 resonance asteroids. Since 
NEAs are often observed at high phase angles this increase in the spectral slopes was interpreted 
as phase reddening 

20101 ) carried out a spectroscopic and spectrophotometric survey of the 



Nathues (200 



Eunomia asteroid family, obtaining spectra in visible (VIS) and NIR wavelengths of 97 of its 
members. The analysis of the spectral slopes of these asteroids revealed an average increase of 
0.067%/100 nm per degree with increasing phase angle in the range of 2° < g < 24°. Apart 
from the increase of the spectral slope he also observed an increase of the 1 (im absorption band 
(de pth) with increasin g phase angle. 

iReddv et al.l (2012j) conducted an extensive study of Asteroid (4) Vesta in order to quantify 
phase angle- induced spectral effects on this asteroid prior to the arrival of the Dawn spacecraft. 
Rotationally resolved NIR spectral observations (0.7-2.5 //m) were obtained for this purpose. 
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They found that in the phase angle range of 0° < g < 25° for every 10° increase in phase angle, 
Vesta's Band I and Band II depths increase 2.35% and 1.5%, respectively, and the BAR value 
increases by 0.30. 

Phase reddening effect has been also detected and quantified from spacecraft observations 
of asteroids. The NEAR Shoemaker spacecraft obtained NIR spectroscopic observations (0.8- 
2.4 fim) of Asteroid (433) Eros at phase angles ranging from 0° to 100°. These observations 
showed the most intense phase reddening for wa yelengths ins i de of the 1.0 //m band, occurring 
at the level of 10% across this phase angle range (jClark et all l2002bl : iBell et all I2002T ) . Similar 
results were observed from NIR reflectance spectra of Asteroid (25143) Itokawa acquired by the 



Hayabusa spacecraft ( Kitazato et al. . 20081 ) 



Laboratory m easurements o f differ ent materials have confirmed the existence of the phase 
reddening effect. Gradie et al. ( 1980l ) investigated phase angle induced effects on the spec- 
trophotometric properties of powdered materials. They obtained reflectance spectra at phase 
angles of 4° < g < 120°, in the wavelength range 0.4 < A < 1.2 fim. All the samples studied by 
them showed a significant reddening as the phase angle increases from 4 to 120°. Furthermore, 
some of the spectra showed variations in the spectral contrast of the absorption bands with 
increasin g phase angle. Mo re re cent laboratory studies relat ed to phase reddening include the 
works of ilYIann et all (|201lh and lShepard and Cloutisl (|201lh . 

Despite the fact that phase reddening has been known for a long time, its effect on the 
analysis of asteroid spectra has not been fully assessed. In the present work we study the 
phase reddening on NEAs and its implications for miner alogical analysis, space weathering 
and taxonomic classification. Due to their proximity, NEAs can exhibit a wide range of phase 
angles, often much higher than those observed for main belt asteroids, which make them the 
logical choice for this study. Our investigation focuses on the analysis of 27 VIS-NIR spectra of 
12 NEAs previously classified as either S-complex or Q-types (Bus-DeMeo), which have been 
observed at different phase angles. In addition to the ground-based observations, laboratory 
spectra of ordinary chondrites are also analyzed and the results compared to those obtained 
from the asteroid spectral data. 



2 Phase reddening from ground-based observations of NEAs 



2.1 The data 

One of the drawbacks of this kind of study is the necessity to obtain multiple observations of the 
same asteroids at different phase angles, which in terms of allocated observa t ion time can be very 
difficu l t to do. For this rea s on we have d e cided to combine data from Reddv ( 20091 ). DeMeo et al. 



d2009l ). lAbell et all (l2007h . lBinzel et al.1 (|200lh and the MIT-UH-IRTF Joint Campaign for NEO 
Spectral Re connaissance ( NEOS R). All of these observations were conducted with the SpeX 
instrument ( Ravner et al. . 20031 ) on the NASA Infrared Telescope Facility (IRTF). The NIR 
spectra (0.7-2.5 //m) were obtained using Spex in its low resolution (R ~ 150) prism mode with 
a 0.8" slit width. Typical observations consist of acquiring frames employing a nodding sequence 
in which the object is alternated between two different slit positions (A and B) following the 
sequence ABBA. Since the object and the sky are observed simultaneously in each exposure, by 
subtracting A from B and B from A, is possible to remove the sky background during the data 
reduction process. In order to minimize the effects of differential atmospheric refraction the 
slit is oriented along the parallactic angle during the observations. To correct for telluric water 
vapor features, and to obtain the relative reflectance, local standard and solar analog stars are 
observed at airmasses as similar as possible as the asteroids. Flat fields and arc line spectra for 
each night are also obtained during the observations. The data reduction is performed using 
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different packages, including IRAF and Spextool (jCushing et al.l. 1200411. F or more detail s abou t 
the observing prot ocols and reductio n of th ese data sets see iBeddvl (|2009l l. lDeMeo et al.l (|2009l l. 
Abell et all (|2007l ) and iBinzel et all (l200lll. The sp e ctra of the NE As at visible wavelen gths 
used in the present work were taken from lBinzel et al.l (j200ll . l2004l l and beMeo et al.1 (l2009l l. All 
the spectra were normalized to unity at 0.55 /im. Observational circumstances for the studied 
objects are presented in Table [TJ 



2.2 Spectral band analysis of NEAs 

As mentioned earlier, increasing the phase angle will produce an increase in the spectral slope 
and changes in the strength of the absorption bands. Spectral slope variations can be inter- 
preted as surface heter ogeneities cau s ed by different factors suc h as metal con tent, particle size 
and space weathering (jGaffev et al.l . Il99.4 Idark et all l2002al : iGaffevl . l20ld l. However, this 
parameter is also known to be very sensitive to other factors. Apart from the viewing geometry, 
atmospheric differential refraction ( Filippenko . 19821 ). air mass differences b etween the standard 



star and the asteroid at the time of the obser vations (jDeMeo et all 120111. incorrect cent ering 



of the object in the slit (jCushing et all I2004T). poor we ather conditions (jBus et all [2002J) and 



the use of different solar analogs ( Sunshine et al. . 20041 ) are often the cause of fluctuations in 
the asteroid's spectral slope. In some cases the origin of these error sources is difficult to deter- 
mine since they are see n as nonsystematic variations in the final asteroid reflectance spectrum 



(|Hardersen et all l200d l. Therefore, in order to better quantify the effect of phase reddening on 



the spectral slope we will leave that part of the study to the laboratory spectra (see section 3), 
and will focus our attention on the analysis of the absorption bands of the asteroid spectra. 

The spectral band parameters for each VIS-NIR spectrum were measur ed using the Spectra l 
Processing Routine (Spe cPR) pr o gram based on the protocols discussed by Cloutis et al. ( 19861 ). 
Gaffev et all (|2002l l. and lGaffev! (|2003l l. Band centers are calculated by dividing out the linear 
continuum and fitting an n-order polynomial over the bottom third of each band. In the case of 
Band II, we defined the unresolved red edge as 2.44 fim. After ratioing the absorption features 
to the straight-line continuum, a subroutine in SpecPR was used to calculate the band areas 
(areas between the linear continuum and the data curve) and band depths (measured from the 
continuum to the band center and given as percentage depths). The band areas were then used 
to calculate the BAR values. Each band parameter was measured ten times using different order 
polynomial fits (typically third and fourth order) and sampling different ranges of points within 
the corresponding intervals. Averages of these measurements were taken as the final values. 
The uncertainties are given by the average 1-a (standard deviation of the mean) calculated from 
the multiple measurements of each band parameter. The olivine-pyroxene abundance ratio of 
the asteroids was estimated using the relationship between oljiol + px) and BAR derived by 
Dunn et al.l (|2010l l from the analysis of 48 ordinary chondrite samples. This linear regression is 
expressed as 



ol/(ol+px) = -0.242 x BAR + 0.728 



(1) 



where the oljipl + px) is expressed as a decimal. The band parameters and ol/(ol + px) 
ratios with their corresponding errors obtained for each asteroid at different phase angles are 
presented in Tables [2] and [3] respectively. 



2.3 Temperature correction to band parameters 

Temperature-induced effects on the spectra of mafic minerals are charac terized by shifting 
the band centers and broadening or narrowing the absorption features ([Singer and Roushl . 
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19851 : iSchade and Waschl . 119991 : iMoroz et all 120001 : Iffinrichs and Lucevl . 120021 ). While sur- 
face temperature variations are considered of minor importance for main- belt asteroid s, they 



may be significant for ob j ects w ith higher eccentric orbits like NEAs (jMoroz et al.l . |2000| ; 



Shestopalov and Golubeva . 2000l ) . Furthermore, the spectra of laboratory samples are com- 
monly obtained at room temperature (~ 300 K). Thus, in order to compare the band parameters 
with those measured for laboratory samples temperature corrections should be applied. 

For each asteroid w e have determined the average surface temperature (T) in the same way 
as (e.g., Burbine et al. . 20091 ). where the temperature of the asteroid is approximated by the 
equation for energy conservation. The calculated temperatures of the NEAs are given in Table 

m 

Since all the asteroids studied in this work are S-complex or Q-types, i.e., olivine-orthopyroxene 
assemblages, we have derived temperature corrections based on the analysis of ordinary chon- 
drites. These meteorites are the most common type of meteorite to fall on Earth ( ~ 85% of all) , 



and they have been general l y linked to either S-com plex or Q-type asteroids (e.g.. iGaffev et al. 



19931 : IVernazza et all 120081 : iNakamura et all l201ll ) 



We have reanalyzed spectra of two H5 ordinary chondrites from iHinrichs and Lucev (|2002l ) 
acquired in the temperature range between 80 and 400 K. Spectral band parameters (band 
centers, band depths and BAR) and their uncertainties were measured using the same meth- 
ods applied t o the asteroid spectra. These data were combined with the results obtained by 
Moroz et al.l (|2000h for a L5 and a LL4 ordinary chondrite acquired at temperatures between 



293 and 80 K. 

We found that Band II is the most affected by temperature variations. A correlation be- 
tween Band II centers and temperature is observed for all samples, i.e., Band II centers shift 
to longer wavelengths as the temperature increases. Therefore, to each data set a linear fit 
was performed and the equations that represent those linear fits were averaged, obtaining the 
following expression: 



BII(T) = 0.0002 x T(K) + 1.87 (2) 

where the Band II center (BII) is given in (im. From this equation we derived a wavelength 
correction for the Band II center 

ABII(nm) = 0.06 - 0.0002 x T(K) (3) 

This correction is derived with respect to room temperature (300 K) and must be added to 
the calculated Band II center of each asteroid. The temperature corrected Band II centers are 
presented in Table [2j For the Band I center we found that the wavelength shift is in general 
very small (~ 0.003 |im), so no temperature correction was derived. 

A decrease in Band II depth with increasing temperature was found for the three groups of 
ordinary chondrites. No obvious trend was observed for the Band I depth. Applying a similar 
procedure to the Band II center, we can estimate the approximate rate of change of the Band 
II depth for a range of temperatures 

BII dep (T) = 20.89 - 0.015 x T(K) (4) 

where the Band II depth (BILi cp ) is given in %. From this equation we derived a temperature 
correction for the Band II depth 

ABII dep (%) = 0.015 x T(K) - 4.5 (5) 

This correction is derived with respect to the room temperature and must be added to the 
calculated Band II depth of each asteroid. 
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An inverse correlation between BAR and temperature was also found for the ordinary chon- 
drites samples. As in the previous cases we obtained an average expression to calculate the rate 
of change of BAR for a range of temperatures 

BAR(T) = 0.83 - 0.00075 x T(K) (6) 
from which a temperature correction with respect to room temperature was derived 

ABAR = 0.00075 x T(K) - 0.23 (7) 

This correction must be added to the calculated BAR values of each asteroid before com- 
paring with those calculated for laboratory samples obtained at room temperature. The tem- 
perature corrected Band II depths and BARs are presented in Table [2j Based on the corrected 
BAR values, the olivine abundance of each asteroid was recalculated using Equation 1. These 
values are included in Table [3l 



T 




Phase Angle (°) 

Figure 1: Band I depth as a function of phase angle for the NEAs. Values derived from multiple 
observations of individual asteroids are connected by lines and have been displayed into two panels for 
clarity. 



2.4 Phase reddening effect on the band parameters 

An inspection of the measured Band I depths presented in Table [2] indicates a change of this 
parameter with phase angle. These values are plotted as a function of phase angle in Figure 
[TJ The black circles represent the measured Band I depths. Values derived from multiple 
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observations of individual asteroids are connected by lines and have been displayed into two 
panels for clarity. A tendency of increasing spectral contrast as the phase angle increases is 
observed for most of the objects in the range of ~ 2° < g < 70°. When higher phase angles are 
included (up to ~ 90°) a slight decrease in Band I depth is observed. However, as we do not have 
sufficient data at such high phase angles, it remains unclear whether this is the actual tendency 
for g > 70°. A similar situation is observed with the temperature corrected Band II depths 
(Fig. [2]). In this case the absorption band seems to increase in the range of ~ 2° < g < 55° 
and then it remains more or less constant for higher phase angles. In some cases we notice a 
deviation from these trends among the measured band depths, i.e., a decrease of band depths 
with increasing phase angle. Although compositional variations cannot be completely ruled 
out, this is more likely due to the fact that some spectra show more scattering than others, 
especially in the Band II. This is probably caused by an incomplete correction of the telluric 
water bands, and due to the decreased response of the detector for wavelengths beyond 2.4 /xm. 
This was seen in the Band I of Asteroids 1036 Ganymed (g = 15.2°) and 4954 Eric (g = 62°) 
and in the Band II of Asteroids 1620 Geographos (g = 38.3°), 1627 Ivar (g = 31.0°) and 1980 
Tezcatlipoca (g = 54.6°). 
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Figure 2: Temperature corrected Band II depth as a function of phase angle for the NEAs. Values 
derived from multiple observations of individual asteroids are connected by lines and have been displayed 
into two panels for clarity. 



In an effort to quantify the phase reddening effect on band depths we performed a linear fit 
to the data of each asteroid. Only the data in the phase angle ranges where the correlations are 
observed were fitted. Those values that deviate from the general tendency were not consider. 
The equations derived from the linear fits were then averaged to obtain a general expression 
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Figure 3: Temperature corrected BAR values as a function of phase angle for the NEAs. Values derived 
from multiple observations of individual asteroids are connected by lines and have been displayed into 
two panels for clarity. 



for each band dep th. Previous works have derived correlations between band parameters and 
phase angle (e.g., Luu and Jewitt . 1990l : Moskovitz et al. . 2O10l ) . However, these correlations 
were obtained from mixed observations of different objects, making it very difficult to disentangle 
the phase reddening from other effects. Our approach, on the other hand, attempts to overcome 
this problem by evaluating the phase reddening on individual objects from which an average 
expression is then derived. We caution that this procedure only provides a rough estimation of 
the effects of phase reddening on band depths, since for each asteroid we have a limited phase 
angle range. The general equations for the band depths are given by 



BI dep (g) = 0.066 x (g) + 17.42 



(8) 



BII dep (g) = 0.093 x (g) + 3.73 (9) 

where the band depths are measured in %. Thus, according to these equations, on average, 
Band I and Band II depths will increase 0.66% and 0.93% respectively, for every 10° increase 
in phase angle in the range of 2° < g < 70° for Band I and 2° < g < 55° for Band II. 

From equations (8) and (9) we can obtain expressions for band depth corrections that are a 
function of the phase angle, 

BI depc = BI dep - 0.066 x (g) (10) 
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Figure 4: Band I center as a function of phase angle for the NEAs. Values derived from multiple 
observations of individual asteroids are connected by lines and have been displayed into two panels for 
clarity. 



BII depc = BII dep - 0.093 x (g) (11) 

where Bl^epc and BIIdepc are the corrected band depths. With these equations we can 
roughly correct the effect of phase reddening in the band depths in the phase angle ranges of 
2° < g < 70° for Band I and 2° < g < 55° for Band II. 

The temperature corrected BAR values presented in Table [2] and plotted in Figure [3] indi- 
cate variations with increasing phase angle, however no obvious trend can be seen. The most 
significant change in BAR is observed for the same objects whose absorption bands show more 
scattering. 

Small variations on band centers with increasing phase angle were found, however they are 
within the uncertainties associated to these band parameters. This can be seen in Figure 2] 
where we plotted Band I center as a function of phase angle for the NEAs. The black circles 
represent the measured Band I centers. Values derived from multiple observations of the same 
asteroids are connected by lines. 
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3 Phase reddening from laboratory measurements of ordinary 
chondrites 



3.1 Data and spectral band analysis 

As it was stated before, ordinary chondrites are considered to be the meteorite analogs of S- 
complex and Q-type asteroids. Therefore, in order to complement our study we have analyzed 
the spectra of a group of ordinary chondrites obtained at a wide range of phase angles. 

Diffuse reflectance spectra were collected at the University of Winnipeg Planetary Spec- 
trophotometer Facility (UWPSF) using an ASD FieldSpec Pro HR spectrometer over the wave- 
length range of 0.35 to 2.5 /im. The three samples that have been analyzed are Dhurmsala 
(LL6, fell 1860), Pavlograd (L6, fell 1826), and Lancon (H6, fell 1897). They were all crushed 
and sieved to a grain size of < 150 [im. The samples were gently poured into aluminum sample 
cups and the edge of a glass slide was drawn across the sample to provide a flat surface for 
the spectral measurements. Reflectance spectra were acquired relative to a 100% Labsphere 
Spectralon disk measured at an incident angle i=13° and emission angle e=0° (13° phase an- 
gle) . Ten sets of measurements were acquired for each sample resulting in three emission angles 
(e=0°, 30°, 60°), five incidence angles (i=0°, 13°, 30°,-30°, 60°) and five different phase angles 
ranging from 13° to 120°. Positive angles (i and e) are measured when the light source and the 
detector are on either side of the normal, while negative incidence angles are measured when 
both light source and detector are on the same side of the normal. For each measurement, a 
total of 250 scans were collected and averaged to improve the signal to noise ratio. 

Spectral band parameters and their uncertainties were measured for each VIS-NIR spectrum 
using the same methods applied to the asteroid spectra. In addition to the band parameters 
we have also measured the spectral slope, which was determined from the fitted continuum 
across Band I, i.e., a straight line tangent to the reflectance peaks from ~ 0.7 to ~ 1.55 /an. 
The uncertainty of the spectral slope is given by the average 1-a, estimated from sampling 
different ranges of points near to the reflectance peaks on either side of the absorption band. 
The olivine- pyroxene abundance ratio of the samples was estimated using Eq.(l). The band 
parameters, spectral slopes and ol/{ol +px) ratios with their corresponding errors obtained for 
each sample are presented in Tabled! 

3.2 Phase reddening effect on the band parameters 

The analysis of the laboratory spectra shows that variations in the band parameters can arise not 
only by changing the phase angle, but also for the same phase angle using different configurations 
of the incidence and emission angle. As can be seen in Table H] phase angles of 30°, 60° and 
90° were obtained using different combinations of i and e, and these different combinations 
produced, in some cases, very different band parameter values. Since we want to quantify the 
effects of phase reddening on the spectral band parameters, for those phase angles where more 
than one combination of i and e was used the average value for each band parameter was taken. 
These average values are presented in Table HI 

Figure [5] shows the reflectance spectra of the LL6 ordinary chondrite at five different phase 
angles. From the bottom to the top the phase angles are 13°, 30°, 60°, 90° and 120°. The 
spectra corresponding to g = 30°, 60° and 90° are the average spectra obtained from different 
combinations of the incidence and emission angles. All the spectra are normalized to unity at 
0.55 ^m. An increase in the spectral slope with increasing phase angle is evident, being more 
significant for phase angles higher than 30°. A similar behavior was observed for the other 
two samples. For phase angles between 13° and 30° a slight reddening is observed in the LL6 
spectra, while a slight blueing (i.e decreasing of the spectral slope with increasing phase angle) 
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is seen in the L6 and H6 spectra. Variations in the strengths of the absorption bands are also 
seen for the three samples. 
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Figure 5: Reflectance spectra of the LL6 chondrite Dhurmsala obtained at five different phase angles. 
From the bottom to the top g=13° (solid line), g=30° (short dashed line), g=60° (long dashed line), 
g=90° (dot-short dashed line) and g=120° (dot-long dashed line). All the spectra are normalized to 
unity at 0.55 /im. 



In Figure [6] we plotted the spectral slopes (panel a), the BAR values (panel b), the Band I 
depths (panel c), the Band II depths (panel d) and the band centers (panels e and f) as functions 
of the phase angle for the LL6 ordinary chondrite. The dashed lines are polynomial fits. 

From the "a" panel we can see that the measured spectral slopes remain more or less constant 
for phase angles between 13° and 30°, and then they increase as the phase angle increases up 
to 120°. A similar trend was observed for the L6 and H6 ordinary chondrites. The largest 
difference in spectral slope between the lowest (13°) and the highest (120°) phase angle was 
found for the LL6 ordinary chondrite, which increased 0.34 /im -1 , followed by the L6 ordinary 
chondrite with an increase of 0.29 //m _1 and the H6 ordinary chondrite with an increase of 0.16 

The Band I depth of the three samples has its maximum spectral contrast near g = 60°. The 
LL6 spectra show a progressive increase of Band I depth from g = 13° to 60° (Fig. [6] panel c), 
however the L6 and H6 spectra show almost no change in Band I depth for phase angles between 
13° and 30°, and then an increase in Band I depth from g = 30° to 60°. The spectra of the 
three samples show a decrease in Band I depth between 60° and 120° phase angle. The largest 
increase in Band I depth between the lowest phase angle (13°) and the phase angle at which the 
Band I reaches its maximum spectral contrast (~ 60° ) corresponds to the LL6 sample, with an 
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Figure 6: LL6 ordinary chondrite: Measured spectral slopes (a), BAR values (b), Band I depths (c), 
Band II depths (d) Band I centers (e) and Band II centers (f) as a function of the phase angle. The error 
bars are smaller than the symbols. The dashed lines represent polynomial fits. 



increase of 2.25%, followed by the L6 and H6 with an increase of 1.36 and 0.53% respectively. 
Of the three samples, the LL6 also shows the most significant decrease of Band I depth between 
the phase angle corresponding to the maximum spectral contrast and the highest phase angle 
(120°), with a decrease of 4.88%, followed by the H6 with a decrease of 4.21% and the L6 with 
a decrease of 3.89%. 

The behavior of the Band II depths seems to be more complex. The LL6 spectra show 
a progressive increase of Band II depths from g = 13° to ~ 45° (where the band reaches its 
maximum spectral contrast) and then it drops between 45° and 90°, becoming more or less 
constant from 90° to 120° phase angle (Fig. [6] panel d). The Band II depths of the L6 spectra 
show almost no change between 13° and 30°. From 30° to 60° Band II depths increase and 
beyond 60° gradually decrease with phase angle increasing to 120°. The Band II depths of the 
H6 spectra show a slight increase from 13° to 30° and then drop from 30° to 120°. The most 
significant increase in Band II depths between the lowest phase angle and the phase angle at 
which the Band II reaches its maximum spectral contrast was observed for the LL6 sample, 
which shows an increase of 1.56%, followed by the L6 with an increase in Band II depth of 
0.56% and the H6 with the lowest increase of 0.16%. The largest difference in Band II depth 
between the phase angle corresponding to the maximum spectral contrast and the largest phase 
angle (120°) was observed for the L6 sample that shows a decrease of 3.11%, followed by the 
H6 with a decrease of 3.02% and the L6 with a decrease of 1.33%. 

Variations in the BAR values are observed among the three samples, being the maximum 
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difference between the lowest and highest phase angle ~0.1, however no systematic trends were 
found. A similar situation occur with the band centers, where the maximum shift (to shorter 
wavelengths) between the lowest and highest phase angle is ~0.01 fim for the Band I center of 
the LL6 sample (Fig. [6] panel e). 



4 Phase reddening and mineralogical analysis 

The analysis of the laboratory samples indicate that the (olivine-rich) LL6 ordinary chondrite is 
the most affected by phase reddening. This sample shows the largest variation in spectral slopes 
and band depths with increasing phase angles. This mat erial-dependence o f the phase reddening 



confi rm the results reported by previous studies (e.g., lAdams and Filicd . 1 19671 : iGradie et al. 



1980h . A comparison between the measurements obtained from the laboratory samples with 



those from the NEAs spectra, reveal that the rate of change in band depths for the ordinary 
chondrites is lower than that found for the asteroids. Differences in grain size, composition and 
packing between the NEAs studied and the meteorite samples could explain this discrepancy. Of 
the three samples, the spectral behavior of the LL6 is the closest to the NEAs. This resemblance 
could be attributed to the fact that half of the asteroids studied are located in the same region 
as the LL ordinary chondrites in the Band I center versus BAR diagram (see Figs. [TOl and [TT|) . 

Because band centers and composition of mafic minerals are related, these band parameters 
are used to derive information about the pyroxene and olivine composition of asteroids surfaces. 
For that purpose, and from the analysis of laboratory samples, different empir ical equations have 
been derived (e.g. , lOaffev et all 120021 : iBurbine et all . 120071 : IPunn et all . l2010h . Our results show 
that there is no significant change in band centers with increasing phase angle for both asteroid 
and meteorite spectra, being the largest shift of ~ 0.01 /jm (for the meteorite spectra), which is 
on the order of the uncertainty associated with this parameter for the asteroid spectra. Based 
on this, it seems to be unlikely that phase reddening could lead to a misinterpretation of the 
minerals composition in asteroids. 

In addition to the band centers, we used the BAR values and Eq. (1) to estimate olivine and 
pyroxene abundances for both NEAs and the laboratory samples. The results obtained (Tables 
OandHD show that the largest variation in the ol/(ol + px) ratio with increasing phase angles 
for NEAs is ~ 0.02 (not taken into account the values obtained from the noisy spectra), and 
for the ordinary chondrites between 0.03 and 0.05. Since these variations are on the order of 
the errors, effects of phase reddening on the estimation of the olivine-pyroxene abundance ratio 
seems to be negligible. 



5 Phase reddening and space weathering 



Space weathering is the term commonly used to refer to any process that modifies the surfaces 
of airless bodies. The effects of space weathering on the spectr a can be seen as reddening 
of the spectral slop es and suppression of the absorption bands ( Pieters et al. . 2000l : Hapkel 



200ll : iGaffevl . l2O10h . The analysis of returned samples from the Moon and Asteroid (25143) 



Itokawa have shown that the cause of the spectral ch anges is the presence of submicroscopic 
metallic iron (SMFe) incorporated into the soil grains (jPieters et all l2000l : iTavlor et all . l200ll 



Noguchi et all . l201lh . This SMFe is produced by condensation of vapors created by microm- 



eteorite impacts and/or deposition of atoms sputter ed off from silicates by solar wind ions 



(ITavlor et all l200ll: iHaoke 



bv lBrunetto and Strazzulla 



2001 



Clark et all l2002ah . Ion irradiation experiments conducted 



(|2005l ) showed that solar wind irradiation can also redden reflectance 



spectra by creation of displacements (the sum of the vacancies and the replacements) caused 



13 



by elastic collisions betw een ions and target nuclei. 

Brunetto etHl (j200fih determined that the effects of space weathering due to ion irradiation 



can be described by an exponential continuum. They computed the ratio between the reflectance 
spectra of irradiated and unirradiated samples and then modeled it with an exponential curve 
given by 



Ratio = W(X) = Kexp(C s /X) (12) 

where A is the wavelength, K is a scale factor and the parameter C5 is a measure of the 
effects of space weathering. They called VF(A) the weathering function. Since both effects, phase 
reddening and space weathering, are manifested in a similar way, we investigated whether the 
red slopes exhibited by spectra obtained at high phase angles could be misinterpreted as space 
weathering. In order to do this we carried out two different experi ments. 

For the first experiment we followed the same procedure used by Brunetto et al. ( 20061 ). with 



the difference that instead of dividing the reflectance spectra of an irradiated by an unirradiated 
sample, we computed the ratio between the spectra obtained at high phase angles and those 
obtained at low phase angles. Since there is no significant change in spectral slope between 
g=13° and 30° we used the spectrum obtained at a phase angle of 30° as reference point. 
Figure [7J shows the ratio plots (solid lines) between the spectra of the LL6 obtained at g=60° 
and g=30° (bottom), at g=90° and g=30° (middle) and at g=120° and g=30° (top pane l ). Th e 
ratio plots were modeled by fitting the same exponential curve used by Brunetto et al. ( 20061 ) 



Ratio = P R (\) = Aexp(C P /X) 



(13) 



where we have just changed the name of the function, we call it the phase reddening function 
(Pr) and the parameters A and Cp are equivalents to the parameters K and C5 respectively, 
with the difference that Cp would be a measure of the effect of the phase reddening. The best fit 
curves obtained using Eq. (13) are represented with dashed lines in Fig. [7J Here we have used 
the LL6 spectra as an example, but the same procedure can be used with the other samples. In 
Fig. [7J we have also included the resultant A and C p values and the coefficient o f determination 
(R 2 ) o f the fitted curves for each ratio plot. Similarly to the results obtained bv lBrunetto et al 



hood ) the ratio of the spectra show some traces of Bands I and II, which can affect the fitting in 
those spectral regions. Nevertheless, the R 2 for the fit ted curves rang e s betw een 0.88 and 0.96 
and, in general, the results are compara ble to those of Brunetto et al. ( 20061 ). After comparing 
the Cp values with the C5 obtained bv iBrunetto et al.1 (J2006J) for different samples, we found 
that the Cp values are similar to the C5 obtained for the H5 ordinary chondrite Epinal, which 
wa s irradiated with Ar ++ (60 keV) at three different ion fluences. 

Brunetto etaD (|2006h found that there is a strong correlation between the C5 parameter 
and the number of displacements per cm 2 (damage parameter). From their experimental data 
they were able to create a damage curve that was fitted by 



C s = aln^d + 1) (14) 

where d is the damage parameter, a = —0.33 ± 0.06 [im and (3 = (1.1 ± 0.5) x 10 -19 cm 2 . 
Using this equation we calculated the number of displacements per cm 2 that would correspond 
to our Cp values. It is important to point out that in this case the d values calculated using 
the C p have no real physical meaning, since the C p parameter quantifies the effect of the phase 
reddening and not the space weathering, however their calculation is useful to illustrate how the 
phase reddening can resemble different degrees of space weathering. This can be seen in Figure 
[8] where we have plotted the Cp (LL6 ordinary chondrite) and C5 (Epinal ordinary chondrite) 
values versus the calculated damage parameter for each sample. The LL6 data are represented 
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Figure 7: Ratio plots (solid lines) between the spectra of the LL6 obtained at g=60° and g=30° (bottom), 
at g=90° and g=30° (middle) and at g=120° and g=30° (top panel). Dashed curves are the best fit 
curves obtained using Eq. (13). For each ratio plot the resultant A, Cp and R 2 values are given. 



with filled symbols while the Epinal values are plotted as open symbols. These results indicate 
that an increase in phase angle in the range of 30° to 120° would be equivalent to irradiate the 
sample with high mass ions (Ar ++ 60 keV) with an ion fluence in the range of ~ 1.3 x 10 15 
to 1.7 x 10 16 ions cm -2 . These ion fluences are equivalent to exposure times of ~ 0.1 x 10 6 to 



1.3 x 10 years at a bout 1 AU from t he Su n (calculated as in lStrazzulla et al.l . l2005l ). Laboratory 



experiments (e.g., Strazzulla et al. . 20051 ) have shown that the timescale for the weathering of 



NEAs surfaces due to ion irradiation is on the order of 10 4 -10 6 years. Thus, observing NEAs 
at high phase angles could produce spectral slopes that resemble those exhibited by weathered 
surfaces. 

For the second experiment we modeled the optic al effec t s of the SMFe on the laboratory 
spectra by using Hapke's radiative transfer model ()Hapkel . Il98ll . Il993l . 1200 li b and then we 
compare d these spectr a with those obtained at different phase angles. 

From lHapke (|200ll l the bidirectional reflectance of a medium of isotropic scatterers ignoring 
the opposition effect is given by 



r(i, 



w 



e,gj 



/'o 



(15) 



47T /i + /U 

where i, e and g are the incidence, emission and phase angle respectively, ^o =cos (i)> //=cos(e), 
wis the single scattering albedo, 7 = (1 — w) 1 / 2 , and H(7, /j) is an analytic approximation to the 
Ambartsumian-Chandrasekhar H functions. Since most reflectances are relative to a standard, 
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Figure 8: The Cp (LL6 ordinary chondrite) and Cs (Epinal ordinary chondrite) parameters as a function 
of the number of displacements per cm 2 (damage parameter). The data corresponding to the Epinal 
meteorite were obtained from iBrunetto et al.l ( 20061 ). The LL6 and the Epinal data are plotted as filled 
and open symbols respectively. For the LL6 each symbol represents the value obtained from the ratiocd 
spectra of different phase angles. The symbols corresponding to the Epinal meteorite represent different 
ions fluences. Approximate exposure times (at 1 AU) corresponding to each ion fluence are also shown. 
The solid line represents the experimental damage curve from Brunetto et al. ( 20061 ). The error bars of 
the LL6 data are on the order of the size of symbols. 



Hapke (l200ll ) represented this relative reflect 



ance as 



rw = (TTWT^, <16) 

from which 7 can be determined, and the single scattering albedo of the sample can be 
calculated as 

w = 1 - 7 2 (17) 

The single scattering albedo can be also written as a function of the properties of the particles 
of the medium by 

w = s e + (i-s e )±^e (is) 

where S e is the Fresnel reflection coefficient for externally incident light, Si is the Fresnel 
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reflection coefficient for internally scattered light (see Hapke . 200 ll ). and is the single-pass 
transmission of the particle. If there is no internal scattering, then 



6 = e- a{D) (19) 

where (D) is the mean ray path length (jHapkd . 119931 ) and a is the absorption coefficient 
given by 



a 



4nnk 
A 



(20) 



where n and k are the real and imaginary part of the refractive index respectively and A is 
the wavelength. From Eqs. (18) and (19) the absorption coefficient can be also calculated as 



1 



a 



(D) 



In 



Si + 



(1-S e )(l-Si 



w 



Sr. 



(21) 



To model the effects of the SMFe, the absorption coefficient of the host material (ah) is 
increased by adding to it the ab sorption coeffi cient of the SMFe (a Fe ). Using the Maxwell- 
Garnett effective medium theory. jHapkd (|200ll ) derived an expression to calculate a Fe , 



OiFe 



\pFe 



(22) 



where / is the mass fraction of the Fe particles, ph is the density of the host material, p Fe is 
the density of iron and z is given by 

= n 3 h n Fe k Fe 

Z ~ (n* Fe - k* Fe + 2n\Y + (2n Fe k Fe y ( ' 

where and n Fe are the real part of the refractive indices of the host material and iron 
respectively, and k Fe is the imaginary part of the refractive index of iron. Thus, if we use 
Eq. (21) to calculate the absorption coefficient of the host material (ah), then the absorption 
coefficient of the material containing SMFe is given by 



a,, 



ah + a Fe 



(24) 



Looking at Eqs. (15) and (18)-(20) we can understand why phase reddening and space 
weathering are manifested in similar ways. Reflectance spectra are controlled by the single 
scattering albedo (w), which is a function of the absorption coefficient (a). The absorption 
coefficient is a parameter that characterizes how deep into a material light of a particular 
wavelength can penetrate before being absorbed. As the phase angle increases, light is less able 
to escape from the medium, meaning that photons are more absorbed, resulting in a decrease 
of reflectance throughout the entire spectrum. However, because the absorption coefficient is 
inversely proportional to the wavelength (the shorter the wavelength the higher the absorption 
coefficient) the reflectance in the blue part of the spectrum will decrease faster than in the red 
part, producing the increase in spectral slope (reddening). Similarly, the addition of the SMFe 
(Eq. 24) will decrease the blue part of the spectrum more than the red part due to the stronger 



absorption at short wavelengths, increasing the spectral slope (jHapkd . 12001 



For our experiment we used the spectrum of the LL6 ordinary chondrite obtained at a phase 
angle of 13° as reference point. Using these data we first calculated ah with Eqs. 16 to 2 1, where 
we have assumed = 1.7, which is a typical value for mafic minerals (jHapkd . l200ll ). Then 
using Eqs. 22 and 23 we estimated a Fe . For th is calculation we have assumed ph = 3.48, which 
is the average grain density for LL chondrites (jBritt and Consolmagn ol 120031 ) and p Fe = 7.87. 
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Figure 9: Reflectance spectra of the LL6 chondrite Dhurmsala after adding different amounts of SMFe 
(dot-short dashed lines) to the spectrum obtained at g=13°. For comparison the spectra of the LL6 
obtained at g=13° (solid lines) and g=60°, 90° and 120° (long dashed lines from the bottom to the top) 
have been included. All the spectra are normalized to unity at 0.55 /im. 



The optical constants of iron, np e and kp e , were taken from Johnson and Christy ( 19741 ). Since 
these values were measured only for A ~ 1.9 /im, we did a polynomial fit to the data in order 
to extrapolate them to 2.5 fim. The ap e was calculated for different mass fractions of Fe and 
then the resulting values were added to the calculated a/j in order to obtain a w . These a w 
values were then inserted into equation (19), which combined with equation (18) allowed us to 
determine the new single scattering albedos. Finally, we calculated 7 from Eq. (17) and using 
Eq. (16) we obtained the "weathered spectra" for the different % of SMFe. In Figure [9] we 
plot the reflectance spectra of the LL6 after adding different amounts of SMFe. We found that 
the spectrum obtained at g = 60° shows a spectral slope comparable to the spectrum obtained 
at g = 13° after adding 0.0015 %SMFe. For phase angles of 90° and 120° the spectral slopes 
are equivalent to adding 0.0048 and 0.0062 %SMFe respectively, to the spectrum obtained at 
13°. If we consider that an amount of ~ 0.02% of SMFe is require d to explain the red 

this means that the 
30% of that 



20041 ). 



slopes of some S-complex asteroids (e.g.. iHapkd . l200ll ; iRivkin et al 
spectral slope for the highest phase angle (g = 120°) would be equivalent to 
SMFe. 
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6 Phase reddening and taxonomic classification 



The taxonomic classification of asteroids is based on shared observational parameters like spec- 
tral slope, color, alb edo and band depth. Most of the curr ent classification systems are based 
on visible data (e.g., Tholen, 1984 ; Bus and Binzell . 2002al lbh . in part because only during the 
last decade sufficient high-quality NIR spectral data became available to extend the classifi- 
cation to the near- infrared wavelengths. Two of th e most common sy stems used to classify 
asteroids using VIS - NIR data are those i ntroduced by Gaffev et al. ( 19931 ) and more recently by 
DeMeo et al. ( 20091 ). Gaffev et al. ( 1993) dey e loped their classification system from the analysis 



of 39 asteroids classified as S-type by iTholenl (|l984h . They divided the S-population into seven 
main compositional subgroups designated S(I)-S(VII). These subgroups range from pure olivine 
through olivine- pyroxene mixtures to pure pyroxene mixtures, and were derived on the basis of 
two band parameters, the Band I center and the BAR. The Bus-DeMeo taxonomy, on the other 
hand, is based on Principal Component Analysis and is comprised of 24 classes that include 
three major complexes (S-, C- and X-complex) and the end members O, Q, R, V, D, K, L, T. 
The S-complex is subdivided into S, Sa, Sq, Sr and Sv. Since both taxonomic systems use band 
depths as one of the primary criteria to classify the objects, multiple observations of the same 
asteroid obtained at different phase angles could lead to ambiguous classifications. In order to 
test the influence of phase reddening on the taxonomic classification, we have applied the two 
classification systems described above to each of the observed asteroids. 

Figure [10] shows the measured Band I center versus BAR for the NEAs. Within the uncer- 
tainties all the NEAs studied are classified either as S (III) or S(IV). There are four cases (1620 
Geographos, 1627 Ivar, 1980 Tezcatlipoca and 1036 Ganymed) that show large variations of the 
BAR from one phase angle to another, but this is likely attributed to noisy spectra rather than 
phase reddening. The rest of the asteroids show variations in Band I centers and BAR that 
are on the order of the error bars. The results obtained from the laboratory samples (Fig. [TTj) 
indicate that the largest variation in the BAR (~ 0.1) is larger than the typical uncertainty 
associated with the BAR values measured from asteroids (represented by the 1-a error bars 
in the upper corner of this figure). This means that in certain cases phase reddening could 
contribute to an ambiguous classification, particularly if the measured band parameters of the 
object are located close to the boundaries that define each class. 

The classification into the Bus-DeMeo system was done using the online taxonomy calcula- 
tor (http://smass.mit.edu/busdemeoclass.html). The class assigned to each spectrum and the 
calculated principal components PCI' and PC2' are given in TableO Only in three cases aster- 
oids were ambiguously classified. Asteroid (1036) Ganymed was classified as S and Sr, (1620) 
Geographos was classifi ed as S and Sq and (4954) Eric was classified as Sw and Sr. The letter 
"w" was introduced bv lDeMeo et al. I (|2009h to indicate an object exhibiting a high slope, but 
does not represent a different class, thus an object designated as "Sw" is an asteroid classified 
as S-type that exhibit a higher s lope than the typica l S-types. This distinction is based on an 



arbitrary cutoff at slope = 0.25 ( DeMeo et al. . 20091 ) . Table [3] shows that the notation of "w" 



was added to many asteroids from one phase angle to another. Certainly the phase reddening 
most play a role on this variation in spectral slope, however as it was stated earlier this is a 
very sensitive parameter that can be affected by other factors. What is important to point out 
is that in the Bus-DeMeo system the spectral slope is removed prior to the classification. For 
those objects which were given two classifications it is likely that the phase reddening was not 
the only contribution, since some of the spectra show more scatter than others. However from 
our analysis we estimate an average variation of about 0.04 and 0.03 for the PCI' and PC2' 
respectively, that could be attributed to phase reddening. These results suggest that phase 
reddening could lead to an ambiguous classification but only if the calculated PC values are 
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close to the lines that separate each class. This can be seen in Figu re 



the PC values in a PCI' versus PC2' diagram from iDeMeo et al.l (|200 



where we have plotted 
). For clarity only the 

calculated PC values for some of the asteroids have been plotted. These values are represented 
by different symbols, which are connected by lines to indicate multiple values for individual 
asteroids. 




Band Area Ratio 



Figure 10: Plot of the Band I center versus BAR for the NEAs. The polygonal region represents the mafic 
silicate components of ordinary chondrites and S(IV) asteroids (jGaffev et al.l . ll993l ). The grey oval region 
represents the mineralogical zone corresponding to the calcic pyroxene-bearing where the S (III) subtypes 

1993h . es the location of the olivinc-orthopyroxene 



are located (IGaffev et al 



mixing line dCloutis et al. . 1986). The top panel shows the measured values for asteroids 1620 Geographos 
(filled circles), 1036 Ganymed (open circles), 1862 Apollo (filled triangles), 1980 Tczcatiploca (crosses), 
11398 (open squares) and 6239 Minos (filled squares). The bottom panel shows the measured values for 
asteroids 4954 Eric (open circles), 66146 (filled squares), 1627 Ivar (crosses), 35107 (open squares), 25143 
Itokawa (filled circles) and 4179 Toutatis (filled triangles). Measured values from multiple observations 
of individual asteroids are connected by lines. The average 1-a error bars are shown in the upper right 
corner. 
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Figure 11: Plot of the Band I center versus BAR for the LL6 (circles), L6 (triangles) and H6 (squares) 
ordinary chondrites. The polygonal region represents the mafic silicate components of ordinary chon- 
drites and S(IV) asteroids ( Gaffev et al. . 19931) . The dashed curve indicates the location of the olivine- 



orthopyroxene mixing line (jCloutis et all 119861 ). The horizontal lines represent the approximate bound- 
aries that separate the three types of ordinary chondrites found by iDunn et al. I (120101) . The error bar in 
the x-axis is of the order of the size of symbols. For comparison purpose we have included the average 1-cr 
error bars (upper right corner) associated with the Band I center and BAR value measured for asteroids. 
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Figure 12: Calculated PC values for some of the NEAs studied represented in a PCI' versus PC2' 
diagram from DeMeo et al. ( 20091) . The letters indicate the different classes within the S-complex plus Q- 
, 0-, R- and V-types. The asteroids whose calculated PC values have been plotted are: (1036) Ganymed 
(open triangles), (1620) Geographos (open circles), (1627) Ivar (filled squares), (4179) Toutatis (filled 
circles), (1862) Apollo (filled triangles), 11398 (open squares) and 66146 (crosses). 
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7 Conclusions 



The analysis of VIS-NIR spectra of 12 NEAs observed at different phase angles has revealed an 
increase of band depths with increasing phase angles in the range of 2° < g < 70° for Band I and 
2° < g < 55° for Band II. Our analysis showed that, on average, Band I and Band II depths will 
increase 0.66% and 0.93% respectively, for every 10° increase in phase angle. From the available 
data we have derived equations that can be used to correct the effects of phase reddening in the 
band depths. Small variations in band centers and BAR values with increasing phase angle were 
also found. Similar trends were observed in the laboratory spectra of three different types of 
ordinary chondrites. In addition, an increase in the spectral slope with increasing phase angles 
was found for the ordinary chondrites. This increase in spectral slope is more significant for 
phase angles higher than 30°. Of the three types of ordinary chondrites the (olivine-rich) LL6 
is the most affected by phase reddening, showing the largest variations in spectral slopes and 
band depths with increasing phase angles. These variations in spectral band parameters seems 
to have no significant impact on the mineralogical analysis though. We have also found that the 
increase in spectral slope caused by phase reddening can mimic the effect of space weathering. 
In particular, an increase in phase angle in the range of 30° to 120° will produce a reddening 
of the reflectance spectra equivalent to exposure times of ~ 0.1 x 10 6 to 1.3 x 10 6 years at 
about 1 AU from the Sun. Furthermore, the increase in spectral slope due to phase reddening 
is comparable to the effects caused by the addition of various amounts of SMFe. These results 
imply that phase reddening should be considered when studying space weathering effects on 
spectral data. Regarding to the taxonomic classification, we found that phase reddening can 
lead to an ambiguous classification, but only if the measured spectral band parameters and/or 
principal components (depending on the taxonomic system) are located close to the boundaries 
that define each class. 
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Table 1: NEAs observational circumstances. The columns in this table are: object number and designa- 
tion, UT date of observation, the phase angle (g) and the heliocentric distance at the time of observation. 



Object 


UT Date 


g (°) 


Heliocentric Distance (AU) 


1036 Ganymed 


09-Mar-05 fe 


5.9 


4.09 


1036 Ganymed 


01-Jun-06 C 


15.2 


3.01 


1620 Geographos 


10-Mar-08 C 


38.3 


1.09 


1620 Geographos 


27-Feb-08 a 


12.0 


1.16 


1620 Geographos 


29-Jan-01 6 


16.7 


1.39 


1627 Ivar 


02-Oct-08 c 


31.0 


1.58 


1627 Ivar 


03-Dec-08 c 


16.0 


1.90 


1862 Apollo 


13-Nov-05 c 


51.7 


1.05 


1862 Apollo 


22-Nov-05 c 


15.4 


1.20 


1980 Tezcatlipoca 


25-Oct-06 c 


54.6 


1.17 


1980 Tezcatlipoca 


05-Jan-07 a 


25.0 


1.48 


4179 Toutatis 


15-Sep-04 C 


27.0 


1.09 


4179 Toutatis 


02-Oct-08 c 


68.0 


1.06 


4179 Toutatis 


13-Aug-08 a 


6.1 


1.45 


4954 Eric 


20-Jul-07 c 


28.3 


1.67 


4954 Eric 


26-Nov-07 a 


62.0 


1.10 


6239 Minos 


26-Jan-04 fe 


47.3 


1.03 


6239 Minos 


06-Sep-10 c 


2.4 


1.20 


11398 


1 n n i nor 

10-Mar-08 


33.0 


1.16 


11398 


28-Feb-08 a 


19.0 


1.21 


25143 Itokawa 


^-Mar-Ol ' 


25.6 


1.07 


25143 Itokawa 


28-Mar-01 e 


69.0 


1.01 


35107 


26-Jul-08 a 


87.0 


1.08 


35107 


27-Dec-02 b 


45.0 


1.20 


66146 


02-Oct-08 c 


70.3 


1.06 


66146 


13-Oct-10 c 


32.1 


1.15 


66146 


06-Sep-10 c 


55.3 


1.16 



a Data from lReddvl (|2009f ) 

b Data from lDeMeo et al.1 (|2009f ) 



c Data from NEOSR, http://smass.mit.edu/minus.html 
d Data fromlAbell et al.1 (|2007l) 



Data from lBinzel et al.1 ( 2001 ) 
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Table 2: NEAs spectral band parameters and their errors. The columns in this table correspond to: 
object number and designation, the phase angle (g), the average surface temperature (T), the Band 

I center (BI±0.01), the Band I depth (BI dcp ±0.3), the Band II center (BII±0.03), the temperature 
corrected Band II center (ABII±0.03), the Band II depth (BIIdep±0.5), the temperature corrected Band 

II depth (ABIIdep±0.5), the band area ratio (BAR±0.04) and the temperature corrected band area ratio 
(ABAR±0.04). 



Object 


g 


T 


BI 


BIdcp 

(%) 


BII 


ABII 


Blldcp 

(%) 


ABIIdcp 


BAR 


ABAR 




(°) 


(K) 


(/jm) 


(/jm) 


(/im) 


(%) 






1036 Ganymed 


5.9 


129.6 


0.909 


18.79 


1.867 


1.901 


8.05 


5.49 


0.901 


0.768 


1036 Ganymed 


15.2 


150.9 


0.917 


18.22 


1.883 


1.913 


9.29 


7.05 


1.101 


0.984 


1620 Geographos 


12.0 


240.4 


0.995 


16.34 


1.970 


1.982 


5.68 


4.79 


0.497 


0.447 


1620 Geographos 


16.7 


219.6 


0.997 


16.77 


1.967 


1.984 


6.56 


5.35 


0.490 


0.425 


1620 Geographos 


38.3 


247.7 


1.001 


18.03 


2.013 


2.023 


4.19 


3.40 


0.278 


0.234 


1627 Ivar 


16.0 


198.9 


0.948 


19.14 


1.955 


1.975 


6.26 


4.74 


0.407 


0.326 


1627 Ivar 


31.0 


218.4 


0.959 


20.28 


1.964 


1.980 


5.04 


3.82 


0.253 


0.186 


1862 Apollo 


15.4 


243.3 


0.986 


25.00 


1.997 


2.008 


4.81 


3.96 


0.188 


0.141 


1862 Apollo 


51.7 


260.0 


0.988 


27.54 


1.989 


1.997 


6.53 


5.93 


0.260 


0.225 


1980 Tezcatlipoca 


25.0 


218.6 


0.947 


16.05 


1.932 


1.949 


6.62 


5.40 


0.526 


0.460 


1980 Tezcatlipoca 


54.6 


246.2 


0.950 


20.45 


1.946 


1.956 


4.37 


3.56 


0.228 


0.183 


4179 Toutatis 


6.1 


222.3 


0.943 


14.91 


1.971 


1.987 


4.89 


3.72 


0.551 


0.487 


4179 Toutatis 


27.0 


256.9 


0.951 


17.44 


1.973 


1.981 


7.29 


6.64 


0.524 


0.487 


4179 Toutatis 


68.0 


260.5 


0.952 


18.04 


1.962 


1.970 


6.90 


6.31 


0.487 


0.452 


4954 Eric 


28.3 


207.1 


0.929 


18.74 


1.925 


1.944 


9.78 


8.39 


1.073 


0.998 


4954 Eric 


62.0 


255.2 


0.923 


16.69 


1.946 


1.955 


8.81 


8.14 


0.917 


0.878 


6239 Minos 


2.4 


244.1 


0.994 


16.56 


2.030 


2.041 


4.03 


3.19 


0.239 


0.192 


6239 Minos 


47.3 


263.2 


0.993 


19.02 


2.010 


2.017 


6.10 


5.55 


0.348 


0.315 


11398 


19.0 


243.3 


0.913 


19.04 


1.935 


1.946 


6.89 


6.04 


0.622 


0.575 


11398 


33.0 


248.3 


0.925 


19.81 


1.945 


1.955 


8.33 


7.55 


0.653 


0.610 


25143 Itokawa 


25.6 


229.2 


0.988 


17.02 


2.017 


2.030 


6.07 


5.01 


0.448 


0.390 


25143 Itokawa 


69.0 


235.2 


0.982 


19.85 


2.017 


2.030 


5.45 


4.48 


0.395 


0.341 


35107 


45.0 


244.1 


0.984 


19.57 


2.019 


2.030 


7.40 


6.56 


0.452 


0.405 


35107 


87.0 


257.6 


0.997 


18.90 


1.996 


2.004 


6.92 


6.28 


0.444 


0.407 


66146 


32.1 


250.1 


0.992 


22.17 


1.983 


1.993 


6.20 


5.45 


0.339 


0.296 


66146 


55.3 


249.1 


0.987 


22.99 


1.984 


1.994 


6.35 


5.59 


0.290 


0.246 


66146 


70.3 


260.2 


0.986 


22.57 


1.992 


2.000 


5.75 


5.15 


0.306 


0.272 
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Table 3: Olivine abundance with its error and taxonomy for the NEAs. The columns in this ta- 
ble are: object number and designation, the phase angle (g), the olivine- pyroxene abundance ratio 
(ol/(ol+px)±0.03), the temperature corrected olivine-pyroxene ratio (Aol/(ol+px)±0.03), taxonomic 
classification (Bus-DcMeo) and principal components PCI' and PC2' for each asteroid. 



Object 


g (°) 


ol/(ol+px) 


Aol/(ol+px) 


Taxonomy 


per 


PC2' 


f036 Ganymed 


5.9 


0.51 


0.54 


S 


0.2400 


0.0416 


1036 Ganymed 


15.2 


0.46 


0.49 


Sr 


0.3957 


0.0381 


1620 Gcographos 


12.0 


0.61 


0.62 


Sq 


-0.1096 


0.0038 


1620 Gcographos 


16.7 


0.61 


0.63 


S 


-0.0169 


0.0206 


1620 Gcographos 


38.3 


0.66 


0.67 


Sqw 


-0.2337 


0.0503 


1627 Ivar 


16.0 


0.63 


0.65 


S 


-0.0151 


0.0517 


1627 Ivar 


31.0 


0.67 


0.68 


Sw 


-0.0591 


0.0689 


1862 Apollo 


15.4 


0.68 


0.69 


Q 


-0.5926 


0.1996 


1862 Apollo 


51.7 


0.66 


0.67 


Q 


-0.6011 


0.2387 


1980 Tezcatlipoca 


25.0 


0.60 


0.62 


Sw 


0.1132 


0.0508 


1980 Tezcatlipoca 


54.6 


0.67 


0.68 


Sw 


-0.0238 


0.0588 


4179 Toutatis 


6.1 


0.59 


0.61 


Sqw 


-0.2825 


0.0972 


4179 Toutatis 


27.0 


0.60 


0.61 


Sq 


-0.3115 


0.1622 


4179 Toutatis 


68.0 


0.61 


0.62 


Sq 


-0.3278 


0.1282 


4954 Eric 


28.3 


0.47 


0.49 


Sw 


0.1552 


0.0433 


4954 Eric 


62.0 


0.51 


0.52 


Sr 


0.4719 


-0.0369 


6239 Minos 


2.4 


0.67 


0.68 


Sq 


-0.2907 


0.0557 


6239 Minos 


47.3 


0.64 


0.65 


Sqw 


-0.3090 


0.1231 


11398 


19.0 


0.58 


0.59 


Sr 


0.0631 


0.0978 


11398 


33.0 


0.57 


0.58 


Sr 


0.1312 


0.0912 


25143 Itokawa 


25.6 


0.62 


0.63 


Sq 


-0.3425 


0.0946 


25143 Itokawa 


69.0 


0.63 


0.65 


Sqw 


-0.3283 


0.1579 


35107 


45.0 


0.62 


0.63 


Sq 


-0.2848 


0.1281 


35107 


87.0 


0.62 


0.63 


Sq 


-0.3077 


0.1182 


66146 


32.1 


0.65 


0.66 


Q 


-0.5262 


0.1994 


66146 


55.3 


0.66 


0.67 


Q 


-0.5737 


0.1824 


66146 


70.3 


0.65 


0.66 


Q 


-0.4978 


0.1938 



26 



Table 4: Spectral band parameters of the ordinary chondrites. The columns in this table correspond 
to: sample type, the incidence angle (i), the emission angle (e), the phase angle (g), the spectral slope 
±0.004, the Band I center (BI±0.003), the Band I depth (BI dcp ±0.1), the Band II center (BII±0.005), 
the Band II depth (BIId cp ±0.2), the band area ratio (BARztO.Ol) and the olivine-pyroxene abundance 
ratio (ol/(ol+px)±0.03). 



Sample 


i 


e 


o 


Slope 


BI 




BII 




BAR 


ol/(ol+px) 




(°) 


(°) 


(°) 


(U1H } 


fum) 

\r J 




(uwl) 

\r J 


(%) 






LL6 


13 





13 


-0.127 


0.984 


36.81 


1.963 


17.27 


0.58 


0.59 


LL6 


30 





30 


-0.094 


0.982 


38.41 


1.959 


18.35 


0.58 


0.59 


LL6 


o 


30 


30 


-0.083 


0.978 


38.47 


1.961 


18.73 


0.61 


0.58 


LL6 


-30 


60 


30 


-0.179 


0.978 


36.69 


1.964 


18.36 


0.63 


0.58 


LL6 


30,0,-30 


0,30,60 


30 


-0.119 


0.979 


37.86 


1.961 


18.48 


0.61 


0.58 


LL6 


60 





60 


0.002 


0.981 


40.28 


1.957 


18.55 


0.55 


0.59 


LL6 





60 


60 


-0.097 


0.969 


37.30 


1.965 


18.89 


0.68 


0.56 


LL6 


30 


30 


60 


-0.043 


0.982 


39.61 


1.960 


18.22 


0.56 


0.59 


LL6 


60,0,30 


0,60,30 


60 


-0.046 


0.977 


39.06 


1.961 


18.55 


0.60 


0.58 


LL6 


60 


30 


90 


0.099 


0.973 


37.91 


1.963 


18.75 


0.62 


0.58 


LL6 


30 


60 


90 


0.184 


0.980 


36.84 


1.961 


16.46 


0.53 


0.60 


LL6 


60,30 


30,60 


90 


0.141 


0.976 


37.38 


1.962 


17.61 


0.57 


0.59 


LL6 


60 


60 


120 


0.212 


0.972 


34.18 


1.959 


17.50 


0.67 


0.57 


L6 


13 


o 


13 


-0.102 


0.949 


42.37 


1.944 


18.70 


0.57 


0.59 


L6 


30 





30 


-0.090 


0.951 


42.71 


1.943 


18.92 


0.56 


0.59 


L6 





30 


30 


-0.089 


0.950 


42.39 


1.945 


19.03 


0.58 


0.59 


L6 


-30 


60 


30 


-0.192 


0.950 


41.58 


1.947 


18.00 


0.57 


0.59 


L6 


30,0,-30 


0,30,60 


30 


-0.124 


0.951 


42.23 


1.945 


18.65 


0.57 


0.59 


L6 


60 


o 


60 


-0.084 


0.954 


43.16 


1.945 


18.21 


0.51 


0.60 


L6 





60 


60 


-0.085 


0.949 


43.36 


1.948 


18.73 


0.59 


0.59 


L6 


30 


30 


60 


0.004 


0.947 


44.68 


1.939 


20.86 


0.62 


0.58 


L6 


60,0,30 


0,60,30 


60 


-0.055 


0.950 


43.73 


1.944 


19.26 


0.57 


0.59 


L6 


60 


30 


90 


0.097 


0.945 


42.53 


1.941 


17.21 


0.51 


0.61 


L6 


30 


60 


90 


0.177 


0.944 


41.69 


1.948 


17.24 


0.56 


0.59 


L6 


60,30 


30,60 


90 


0.137 


0.944 


42.11 


1.945 


17.23 


0.54 


0.60 


L6 


60 


60 


120 


0.191 


0.944 


39.84 


1.943 


16.15 


0.51 


0.61 


H6 


13 





13 


-0.141 


0.934 


35.45 


1.930 


18.41 


0.88 


0.52 


H6 


30 





30 


-0.151 


0.935 


36.20 


1.929 


19.12 


0.89 


0.51 


H6 





30 


30 


-0.115 


0.933 


35.33 


1.934 


18.33 


0.88 


0.52 


H6 


-30 


60 


30 


-0.177 


0.931 


34.80 


1.931 


18.26 


0.93 


0.50 


H6 


30,0,-30 


0,30,60 


30 


-0.148 


0.933 


35.44 


1.932 


18.57 


0.90 


0.51 


H6 


60 





60 


-0.149 


0.935 


36.48 


1.926 


18.91 


0.83 


0.53 


H6 





60 


60 


-0.096 


0.929 


34.54 


1.935 


16.93 


0.84 


0.52 


H6 


30 


30 


60 


-0.096 


0.933 


36.93 


1.931 


17.85 


0.76 


0.55 


H6 


60,0,30 


0,60,30 


60 


-0.114 


0.933 


35.98 


1.931 


17.90 


0.81 


0.53 


H6 


60 


30 


90 


-0.058 


0.932 


34.77 


1.928 


17.65 


0.83 


0.53 


H6 


30 


60 


90 


-0.012 


0.926 


31.69 


1.934 


14.77 


0.84 


0.53 


H6 


60,30 


30,60 


90 


-0.035 


0.929 


33.23 


1.931 


16.21 


0.83 


0.53 


H6 


60 


60 


120 


0.020 


0.932 


31.77 


1.932 


15.55 


0.80 


0.54 
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A Asteroid spectra 



This appendix contains the spectra of all NEAs analyzed in this study. The numerical designa- 
tion and the date of the observation (YYMMDD) for each asteroid are given. 




Figure 13: Reflectance spectra of NEAs analyzed in this study. 
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Figure 14: Reflectance spectra of NEAs analyzed in this study. 
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Figure 15: Reflectance spectra of NEAs analyzed in this study. 
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